AD-A098  422 
UNCLASSIFIED 


ROCKWELL  INTERNATIONAL  THOUSAND  OAKS  CA  SCIENCE  CENTER  F/S  11/2 
RESEARCH  OF  MICROSTRUCTURALLY  DEVELOPED  TOUGHENING  MECHANISMS  I—ETCtU> 
MA?  6±  £  F  LANGE  N00014-77-C-0441 

SC5il7«9TR  ^ 


ADA098  422 


A 

TRANSFORMATION  TOUGHENING 
PART  2:  CONTRIBUTION  TO  FRACTURE  TOUGHNESS 

F.F.  Lange 

Structural  Ceramics  Group 
Rockwell  International  Science  Center 
Thousand  Oaks,  California  91360 


rtf** 


.  t-i 


ABSTRACT 


Two  approaches  are  taken  to  determine  the  contribution  of  a  stress- 
induced  phase  transformation  to  the  fracture  toughness  of  a  brittle  material. 
Both  approaches  result  In  the  expression: 
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where  K0  Is  the  critical  stress  Intensity  for  the  material  without  the  trans¬ 
formation  phenomenon,  (UGC|  -  Usef)  Is  the  work  done  per  unit  volume  by  the 
stress  field  to  Induce  the  transformation,  Ec  and  vc  are  the  elastic  properties, 
Vj  Is  the  volume  fraction  of  retained,  high  temperature  phase  and  R  Is  the  size 
of  the  transformation  zone  associated  with  the  crack.  It  Is  assumed  that  only 
those  Inclusions  (or  grains)  close  to  the  crack's  free  surface  will  contribute 
to  the  fracture  toughness;  thus  R  *  the  Inclusion  size.  The  chemical  free 
energy  change  associated  with  the  transformation  (|aGc|)  will  govern  the 
temperature  and  alloying  dependence  of  the  fracture  toughness. 
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where  K.  Is  the  critical  stress  Intensity  for  the  material  without  the  transVl 
formation  phenomenon,  (I  AGC I  -  Usef)  Is  the  work  done  per  unit  volume  by  the^S 
stress  field  to  Induce  the  transformation,  Ec  and  vc  are  the  elastic  proper-  '  „ 
ties,  Is  the  volume  fraction  of  retained,  high  temperature  phase  and  R  is 
the  size  of  the  transformation  zone  associated  with  the  crack.  It  Is  assumed 
that  only  those  inclusions  (or  grains)  close  to  the  crack's  free  surface  will 
contribute  to  the  fracture  toughness;  thus  R  91  the  inclusion  size.  The 
chemical  free  energy  change  associated  with  the  transformation  (I  AGC I  will 
govern  the  temperature  and  alloying  dependence  of  the  fracture  toughness. 
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1.  Introduction 

The  first  part^  of  this  series  described  the  thermodynamics  of  the 
constrained  phase  transformation  with  particular  emphasis  on  the  effect  of  the 
Inclusion  size.  In  this  part,  two  different  approaches,  viz.  Griffith's  and 
Irwin's,  will  be  used  to  determine  the  contribution  of  a  stress- Induced  phase 
transformation  to  the  fracture  toughness  of  a  brittle  material.  As  In  Part  1, 
special  reference  will  be  made  to  the  Zr02  ( tetragonal ) ^ Zr02  (nonocllnlc) 
transformation. 
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2.  Griffith’s  Approach 

In  this  approach, ^  the  total  energy  of  a  system,  defined  as  a  cracked 
body  and  the  applied  load.  Is  determined  as  a  function  of  the  crack  area.  The 
critical  condition  for  crack  extension,  as  first  defined  by  Griffith,  corre¬ 
sponds  to  the  maximum  In  the  total  energy  vs  crack  area  function.  To  carry  out 
this  analysis,  we  first  determine  and  sum  the  different  contributions  to  the 
total  energy  that  change  as  a  function  of  crack  area,  viz.  the  body's  strain  and 
surface  energies  and  the  load's  potential  energy.  The  condition  for  crack 
extension  Is  then  found  by  determining  the  maximum  In  the  total  energy  vs  crack 
area  function.  For  the  needs  of  the  present  analysis ,  we  can  simplify  the 
analysis  by  choosing  a  crack  and  loading  system  In  which  all  of  the  required 
functions,  except  the  one  associated  with  the  stress-induced  transformation,  are 
already  known.  The  penny- shaped  crack  under  an  applied  tensile  load  first 
analyzed  by  Sack^  was  chosen  fdr  this  analysis. 
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Figure  1  Illustrates  a  section  of  the  cracked  body  under  tensile  load  ' 

i 

j 

which  contains  a  uniform  dispersion  of  untransformed  Inclusions  of  volume  frac- 

tlon  Vf.  It  Is  assumed  that  the  stress  field  associated  with  the  crack  front  j 

has  caused  Inclusions  to  transform.  As  the  crack  extends  and  the  stresses  j 

within  the  previously  transformed  zone  decrease.  Inclusions  that  have  lost  some 

constraint  by  being  either  traversed  by  the  crack  or  In  close  proximity  to  the 

new  fracture  surfaces  will  remain  In  their  transformed  state.  This  process 

leads  to  a  transformation  zone  which  surrounds  the  crack,  as  shown  In  Fig.  1. 

Only  those  Inclusions  which  remain  transformed  will  contribute  to  the  non- 
recoverable  work  done  by  the  loading  system  to  stress-lnduce  the  transforma¬ 
tion.  Since  the  Inclusion  must  lose  constraint  to  remain  In  their  transformed 
state  once  the  crack's  stress  field  moves,  the  size  of  the  transformed  zone  (R) 

Is  approximately  equal  to  the  Inclusion  size,  viz.  R  *  0. 

Following  Sack's  solution  to  this  same  problem  without  the  transforma¬ 
tion  phenomena,  the  Increase  In  free  energy  of  the  body  due  to  new  crack  sur¬ 
faces  Is 

vA.  .  in 


where  G0  Is  the  critical  strain  energy  release  rate  associated  with  the  forma¬ 
tion  of  new  surface.  The  Increase  In  strain  energy  due  to  crack  extension  Is 
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Fig.  1  Section  of  a  penny-shaped  crack  of  radius  c  within  a  stressed  composite 
containing  Inclusions.  Hatch  Inclusions  have  been  transformed  during 
crack  extension. 
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where  aa  Is  the  applied  tensile  stress  and  Ec  and  vc  are  Young's  Modulus  and 
Poisson's  ratio  of  the  composite  material,  respectively.  The  work  done  by  the 
loading  system,  l.e.,  the  decrease  In  the  potential  energy  of  the  load.  Is 
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An  additional  term  arises  when  the  transformation  phenomena  Is  Included 
In  the  fracture  process.  Recognizing  that  the  volume  of  the  transformed  zone  Is 
2vR(c  +R)2  3  2wRc2,  this  additional  term  Is  the  work  done  by  the  loading  system 
to  form  the  transformed  zone: 


W2  ■  -2irRc2HVi 

where  W  Is  the  work  per  unit  volume  of  transformed  material  to  induce  the  trans¬ 
formation.  The  minimus  value  of  W  Is  determined  by 


'  -«GC  ♦  '  "  '  0  >  141 

where  &GC  Is  the  chemical  free  energy  change  for  the  reaction  ZrOj  (tetragonal) 
♦  Zr02  (monocllnlc),  MJ$e  Is  the  change  In  strain  energy  associated  with  the 
transformation  and  (1  -  f)  Is  the  loss  of  strain  energy  due  to  the  loss  of 
constraint  Imposed  on  the  Inclusions  during  crack  extension.  Thus, 

H  -  -  AG  c+  au$ef  (5) 
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where  KQ  *  [(G0EC)/(1  -  Is  the  critical  stress  Intensity  factor  of  the 

material  without  the  transformation  phenomenon. 

*The  absolute  brackets  are  used  to  Indicate  that  the  sign  of  has  already 
seen  defined  as  negative  In  Eq.  (4)  for  the  temperature  range  of  Interest. 
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3.  Irwin's  Approach 

The  same  subject  can  be  viewed  with  Irwin's  approach/*)  where  one 
Imposes  a  force  field  to  a  stressed  crack  tip  and  calculates  the  work  required 
to  close  It  by  a  unit  length.  Irwin  showed  that  the  work  per  unit  length  of 
crack  closure  Is  equivalent  to  the  net  work  dissipated  per  unit  length  of  crack 
extension  3{Use  -  Wjl/ac  as  calculated  through  the  Griffith  approach.  Fracture 
will  take  place  when  3(U$e  -  Wj)/3C  >  Gc,  termed  the  critical  strain  energy 
release  rate.  The  major  difference  between  the  two  approaches  Is  that  In 
calculating  the  net  work  dissipated,  Irwin  only  needs  to  consider  the  stress 
field  In  the  vicinity  of  the  crack  tip,  whereas  the  Griffith  approach  requires 
complete  knowledge  of  the  stress  state  in  the  system. 

To  apply  the  Irwin  approach,  let  us  consider  the  unit  of  crack  to  be 
closed  to  have  traversed  a  transformed  Inclusion  as  shown  In  Fig.  2(a).  The  work 
to  dose  this  unit  of  crack  can  be  broken  Into  two  parts,  one  concerned  with  the 
transformation! AWt)  and  one  concerned  with  crack  closure  ( aWc) .  The  first  force 
field  we  apply  would  revert  the  fractured  Inclusion  back  to  its  untransformed 
state  as  shown  In  Fig.  2(b)  The  work  performed  by  this  first  force  field  per 
unit  volume  of  transformed  material  Is 


W  »  |A  Gc|  -  AU$ef  (10) 

where  IaGc|  Is  the  change  In  chemical  free  energy  for  reverting  the  Zr02  Inclu¬ 
sion  from  Its  monoclinic  structure  to  Its  tetragonal  structure.  AU^  Is  the 
strain  energy  associated  with  the  transformation,  and  (1  -  f)  Is  that  portion  of 
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the  strain  energy  relieved  during  fracture.  The  total  work  done  on  all 
Inclusions  within  the  volume  2Rac  per  unit  crack  length  Is 

AWt  =  2RACV.W  *  2RVi(|AGC|  -  AU^f)  AC  (11) 


Once  the  inclusions  have  been  reverted  to  their  untransformed  state, 
the  strain  energy  associated  with  the  Inclusions  disappear  and  the  crack  now 
looks  like  any  ordinary  crack  In  a  two  phase  material.  At  this  point,  the 
second  force  field  can  be  applied,  as  defined  by  Irwin,  to  close  the  crack  by 
the  unit  length  ac  as  shown  In  Fig.  2(c).  The  work  performed  In  this  operation 
is 


AWC  =  Gqac  ,  (12) 

where  G0  Is  the  critical  strain  energy  release  rate  for  the  composite  material 
without  the  transformation  phenomenon. 

The  total  work  for  crack  closure  per  unit  crack  length  which  also 
reverts  the  Inclusions  to  their  Initially  untransformed  state  Is 


AW 

AC 


AW£+  AW, 
AC  AC 


G0  +  2RV j( I aGC I  -  AUsef) 


(13) 


Thus,  the  contribution  of  the  stress- induced  transformation  to  the  critical 
strain  energy  release  rate  of  the  composite  can  be  expressed  as 
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Gc  =  G0  +  2RV1( I AGC1  -  AUsgf) 


(14) 


or  expressed  as  the  critical  stress  intensity  factor 


EcGc  1 

1/2 

a 

It  ,  2«»,ECII«5CI  -  *1 „« 

(1  -  v;) 

0  a  -  T) 

c 

c  J 
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(15) 


4.  Discussion 

As  expected,  both  approaches  have  led  to  the  same  expression  for  the 
critical  stress  Intensity  factor  of  a  material  containing  inclusions  that  can 
undergo  a  stress- induced  transformation.  The  expressions  show  that  the  contribu¬ 
tion  of  the  stress  Induced  transformation  can  be  maximized  by  maximizing  a)  the 
volume  fraction  (V^)  of  inclusions  fabricated  In  their  untransformed  state,  b) 
the  elastic  modulus  (Ec)  of  the  composite,  c)  the  factor  (|aGc|  -  AUsef)  and,  d) 
the  size  of  the  transformation  zone  (R)  associated  with  the  propagating  crack 
front.  Each  of  these  factors  will  be  discussed  In  the  following  paragraphs. 

Maximizing  the  volume  fraction  would  result  in  a  single  phase,  poly¬ 
crystalline  material,  e.g.,  tetragonal  Zr02.  Here,  the  Inclusions  could  be 
defined  as  Individual  grains,  surrounded  by  neighboring  grains  of  different 
mlsorlentatlons  which  define  the  matrix.  Each  neighboring  grain  constrains  one 
another  from  undergoing  a  stress-free  transformation  by  their  anisotropic 
transformation  strains.  Thus,  from  a  conceptlonal  viewpoint,  a  single  phase, 
polycrystalline  material  can  be  treated  In  the  same  manner  as  described  above. 
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The  elastic  modulus  of  the  composite  can  be  Increased  by  choosing  a 
chemically  compatible  second  phase  with  a  higher  elastic  modulus.  For  the  case 
of  ZrOg,  AI2O3,  with  a  modulus  approximately  twice  that  of  Zr02,  would  be  a  de¬ 
sirable  choice.  But  adding  a  second  phase  to  increase  the  modulus  would  at  the 
same  time  decrease  the  volume  fraction  (Vj)  of  the  toughening  agent.  One  would, 
therefore,  be  concerned  with  the  product  of  VjEc  (Eq.  (15))  In  optimizing  Kc. 

If  we  assume  that  the  composite  modulus  is  governed  by  the  rule  of  mixtures, 
viz.  E£  =>  EjVj  +  Em(l  -  V j) ,  then  the  product  would  have  the  following  form: 

V1Ec  =  V1E1[M  '  V1(M  ‘  11]  •  (l6) 


where  M  =  Em/Ej,  the  modular  ratio  of  the  matrix  (m)  and  the  Inclusion  (1). 
Differential ■*;  Eq. ( 16)  with  respect  to  Vj, 


3(V .E  ) 

■  e^Cm  -  zv^m  -  1)] 


(17) 


shows  that  the  maximum  In  the  product  relation  occurs  at  Vj  *  1  when  M  <  2, 
l.e.,  the  greatest  toughness,  with  other  factors  constant,  should  be  obtained 
for  a  single  phase  material.  If  M  >  2,  the  toughness  could  be  optimized  when 
Vj  <  1.  On  the  other  hand.  If  the  objective  Is  to  toughen  a  matrix  phase 
(e.g.,  toughening  A^Oj  with  Zr02),  Eq.  (16)  shows  that  more  toughening  Is 
obtained  for  a  given  volume  fraction  the  greater  the  modular  ratio.  That  Is, 
small  volume  fractions  (e.g.,  <  0.3)  of  the  toughening  agent  will  produce 

greater  results  the  greater  the  modulus  of  the  matrix  material. 
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The  dependence  of  fracture  toughness  on  temperatures  and  alloy  content 
Mill  be  governed  by  the  factor  (l&Gcl  -  AUsef).  This  Is  because,  relative  to 
other  factors,  the  chemical  free  energy  change! AGC)  exhibits  the  greatest  change 
Mlth  temperature  and  alloying  content.  For  the  Zr02(t)  ♦  Zr02(m)  reaction, 

I aGc |  decreases  with  Increasing  temperature  and  alloying  (e.g.,  Y2O3,  Ce02,  etc) 
content.  Thus,  for  this  transformation,  the  fracture  toughness  Is  expected  to 
decrease  with  Increasing  temperature  as  the  factor  (|aGc|  -  AUsef)  decreases  to 
zero.  Likewise,  Kc  will  decrease  as  the  alloy  content  In  Zr02(t)  Is  Increased. 
That  Is,  fracture  toughness  will  be  optimized  at  the  lower  temperatures  and  for 
the  least  alloy  content.  The  temperature  where  the  contribution  of  the  stress- 
induced  toughness  disappears  will  depend  on  the  magnitude  of  AUsef.  Phenomena 
that  help  relieve  strain  energy  during  the  fracture,  e.g.,  twinning,  will 
decrease  the  value  of  f  and  thus  Increase  the  temperature  where  the  contribution 
to  toughness  disappears. 

The  major  assumption  used  In  the  model  to  derive  the  Kc  expressions  was 
that  the  size  of  the  transformation  zone  (R)  was  determined  by  the  close  prox¬ 
imity  of  the  Inclusions  to  the  free  surface  formed  during  fracture.  That  is. 
Inclusions  transformed  by  the  stress  field  would  only  remain  In  their  trans¬ 
formed  state,  once  the  crack's  stress  field  passes.  If  much  of  their  constraint 
was  lost  during  crack  extension.  This  assumption  leads  to  the  hypothesis  that 
the  zone  size  would  be  directly  related  to  the  Inclusion  size  (0),  viz.  R  ■  0. 

It  Is,  therefore,  hypothesized  that  Kc  will  Increase  with  Increasing  Inclusion 
size. 
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It  can  also  be  argued  that  If  Inclusions  remote  to  the  crack's  surface 
were  to  remain  In  their  transformed  state,  their  residual  strain  energy  would  be 
greater  relative  to  those  adjacent  to  the  crack  surface.  Thus,  the  work  loss  to 
the  fracture  process  for  remote  Inclusions  would  be  less  than  for  Inclusions 
adjacent  to  the  crack  surfaces.  That  Is,  adjacent  Inclusions  would  contribute 
more  to  the  fracture  toughness  than  remote  Inclusions. 
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